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ABSTRACT 

Natural  convection  from  a  single  column  of  eight  in-line,  rectangular,  flush 
heat  sources  was  examined.  A  vertical  channel  was  formed  by  placing  a  movable 
shroud  parallel  to  the  test  surface.  The  experimental  program  consisted  of  temper¬ 
ature  measurements  and  flow  visualizations.  For  a  range  of  component  power  levels 
and  channel  spacings,  surface  temperatures  of  each  heat  source  were  measured  at 
five  locations  using  embedded  thermocouples.  The  resulting  steady  state  data  was 
compared  with  corresponding  measurements  for  an  identical  geometric  arrangement 
of  protruding  heat  sources.  The  flow  visualization  was  carried  out  using  a  laser 
generated  plane  of  light. 
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Natural  convection  from  a  single  column  of  eight  in-line,  rectangular,  flush 
heat  sources  was  examined.  A  vertical  channel  was  formed  by  placing  a  movable 
shroud  parallel  to  the  test  surface.  The  experimental  program  consisted  of  temper¬ 
ature  measurements  and  flow  visualizations.  For  a  range  of  component  power  levels 
and  channel  spacings,  surface  temperatures  of  each  heat  source  were  measured  at 
five  locations  using  embedded  thermocouples.  The  resulting  steady  state  data  was 
compared  with  corresponding  measurements  for  an  identical  geometric  arrangement 
of  protruding  heat  sources.  The  flow  visualization  was  carried  out  using  a  laser 
generated  plane  of  light. 
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I.  INTRODUCTION 


A.  STATEMENT  OF  PROBLEM 

With  each  successive  increase  in  the  scale  of  integration  of  micro-electronic 
circuits,  there  has  been  a  need  to  improve  cooling  technologies.  Currently,  chip  heat 
fluxes  are  approaching  50  W/cm 2,  [Ref.  1].  Heat  dissipation  may  be  improved  by 
using  higher  velocities  of  air,  chilled  water  heat  exchangers,  or  altering  the  system 
design  to  reduce  its  thermal  resistance.  These  techniques  are  becoming  less  satis¬ 
factory  and  the  limits  of  air  and  indirect  liquid  cooling  have  nearly  been  reached.  If 
projected  heat  fluxes  of  200  W/cm2,  [Ref.  1],  are  to  be  achieved,  then  direct  liquid 
immersion  becomes  an  attractive  option,  [Ref.  1,  6]. 

The  performance  and  reliability  of  micro-electronic  devices  are  very  sensitive 
to  junction  temperatures.  The  failure  rate  of  chips  increase  exponentially  as  the 
junction  temperature  increases.  The  long  term  reliability  of  circuit  chips  can  be 
improved  by  50%  with  every  10°C  decrease  in  temperature  that  the  device  is  main¬ 
tained  below  its  maximum  operating  temperature,  [Ref.  2,  3,  4]. 

For  a  complete  thermal  analysis  of  any  electronic  package  design,  all  three 
modes  of  heat  transfer  must  be  considered.  Radiation  effects  are  significant  primarily 
in  systems  involving  air  cooling.  Effects  of  conduction,  however,  must  be  evaluated 
in  each  cooling  application.  These  may  also  be  substantial  during  transients. 

The  convection  heat  transfer  rate  from  a  surface  at  temperature  T,  with  ex¬ 
posed  surface  area  A  to  its  surrounding  fluid  at  temperature  Tamf,  is  given  by  New¬ 
ton’s  Law  of  Cooling: 


q  =  hA(T,  -  Tamb ) 


(1.1) 


This  is  the  definition  of  h,  the  convective  heat  transfer  coefficient  and  not  a  phe¬ 
nomenological  law. 

Most  available  correlations  for  the  convective  heat  transfer  coefficient  are  usu¬ 
ally  valid  for  large  heat  sources.  These  are  often  good  for  determining  the  overall 
heat  transfer  from  a  single  circuit  board.  Individual  devices  mounted  on  the  sub¬ 
strate  act  as  discrete  heat  sources  and  typically  have  higher  average  temperatures 
than  those  predicted  by  such  text  book  correlations  [Ref.  6]. 

Convection  may  be  either  forced  or  natural.  Forced  convection  relies  on  a 
prime  mover  to  displace  the  fluid  and  cause  it  to  flow  over  the  electronic  package. 
This  method  of  convection  yields  very  high  heat  transfer  coefficients  and  may  become 
one  of  the  primary  means  of  cooling  electronic  packages  as  component  heat  fluxes 
continue  to  rise.  Natural  convection  or  buoyancy  induced  flow,  results  from  body 
forces  acting  on  density  differences  within  the  fluid.  In  many  applications,  these 
density  gradients  result  due  to  the  presence  of  temperature  variations.  The  heat 
transfer  coefficients  are  smaller  than  those  associated  with  forced  convection  since 
the  flow  velocities  are  generally  much  smaller  for  natural  convection.  However,  the 
resulting  design  of  the  package  can  be  considerably  simpler  and  more  reliable  than 
in  other  cooling  arrangements. 

B.  OBJECTIVES 

This  experimental  study  was  undertaken  to  determine  the  hydrodynamic  and 
thermal  characteristics  of  buoyancy  induced  flow  resulting  from  several  flush  mounted 


2 


discrete  heat  sources  on  a  vertical  surface.  The  results  from  this  study  will  be  com¬ 
pared  to  those  in  Reference  7  where  an  array  of  rectangular  protruding  heat  sources 
was  used. 

A  single  column  of  discrete  heat  sources  was  chosen  since  it  provided  a  simple 
geometric  arrangement  simulating  some  of  the  modern  electronic  packages  such  as 
flat  packs. 

Specific  objectives  of  this  study  were: 

•  To  visualize  the  steady-state  natural  convection  flow  within  a  channel  over  a 
range  of  component  power  levels  and  channel  widths. 

•  To  measure  component  temperatures  at  each  power  level  and  channel  width 
and  develop  heat  transfer  relationships  for  this  geometry. 

•  To  measure  the  temperature  distribution  within  the  fluid  adjacent  to  the 
heated  surface  in  the  absence  of  a  shrouding  wall. 

This  study  advances  the  work  described  in  References  5  and  7  and  is  part  of 
an  ongoing  effort  to  study  natural  convection  liquid  immersion  cooling  of  electronic 
equipment. 
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II.  EXPERIMENTAL 

A.  GENERAL  DESIGN  CONSIDERATION 

Equipment  design  and  construction  details  are  provided  in  Reference  5.  Only 
system  changes  and  a  brief  overview  is  presented  here. 

As  seen  in  Figure  2.1,  the  test  assembly  consists  of  eight  flush  mounted  foil 
heaters.  Each  heater  (Figure  2.2)  is  a  network  of  Inconel  600  conductor  mounted  on 
a  Kapton  backing.  The  heaters  measured  23.9  mm  in  width  by  7.8  mm  in  height. 
They  formed  a  single  column  on  a  plexiglass  substrate.  The  substrate  measured 
25.4  mm  square  and  was  12.7  mm  thick.  The  column  of  heaters  was  centered  on  the 
substrate  and  each  heater  was  spaced  25.4  mm  center  to  center  from  an  adjacent 
heater.  Another  piece  of  plexiglass  identical  to  the  substrate  without  the  heaters 
was  used  as  an  adjustable  shroud  wall  to  form  the  channel.  This  shroud  was  used 
to  simulate  the  back  side  of  a  neighboring  circuit  board. 

Five  copper-constantan  thermocouples  of  0.127  mm  diameter  wire  were  used 
for  temperature  measurements  on  each  heater.  The  thermocouples  were  embedded 
within  6  mm  deep  slots  in  the  substrate.  For  each  heater,  one  thermocouple  bead 
was  located  at  the  foil  center  and  the  remaining  four  were  at  the  centers  of  each  edge 
of  the  heater.  All  wires  were  routed  along  6  mm  deep  grooves  cut  below  the  heater 
placement  slots.  All  wiring  came  out  through  the  back  of  the  test  surface  through 
drilled  holes  in  the  substrate.  After  all  heaters  and  wires  were  in  place,  these  holes 
were  sealed  with  a  low  conductivity  adhesive.  Omega  Bond  101,  a  high  conductivity 
adhesive,  was  used  to  bond  the  thermocouples  in  place  and  to  electrically  insulate 
the  thermocouple  beads  from  the  heater  foils.  A  sixth  thermocouple  for  each  heater 
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Figure  2.1:  Assembled  Test  Surface 

was  centered  behind  the  heater  on  the  back  side  of  the  plexiglass  substrate.  This 
was  used  to  estimate  conduction  losses  through  the  substrate. 

A  0.127  mm  diameter,  copper-constantan  thermocouple  probe  (Figure  2.3)  was 
constructed.  This  probe  was  attached  to  a  three  dimensional  traverse  and  used  to 
measure  the  temperature  distribution  within  the  fluid  adjacent  to  the  test  surface. 
The  foil  heaters  were  powered  using  a  regulated  D.C.  power  supply  as  described  in 
Reference  5. 

The  test  surface  was  immersed  in  an  open,  one  cubic  meter  tank  constructed 
from  plate  glass  and  insulated  on  all  sides.  For  flow  visualization,  particles  of  pliolite 
were  suspended  in  the  de-ionized  water.  Pliolite  is  an  inert  substance  with  a  specific 
gravity  of  0.93.  These  particles  were  illuminated  by  aiming  an  eight  milliwatt  helium 
neon  laser  through  a  cylindrical  lens  to  split  the  light  beam  into  a  plane  of  light. 
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Figure  2.2:  Foil  Heater 


Figure  2.3:  Traversing  Thermocouple  Probe 

By  using  a  motor  driven  camera  and  an  intervalometer,  extended  time  exposure 
photographs  of  the  flow  patterns  within  the  illuminated  plane  were  obtained. 

B.  EXPERIMENTAL  PROCEDURE 

Details  of  data  acquisition  techniques  are  provided  in  Reference  5,  hereafter 
referred  to  as  the  initial  study.  The  follow-on  study  by  Willson  [Reference  7]  will  be 
referred  to  as  the  previous  study. 

Power  settings  of  0.2,  0.5,  1.0,  1.2,  1.5,  and  2.0  watts  were  used  for  each 
channel  width.  The  shroud  was  adjusted  for  channel  spacings  of  1.5  mm,  3.0  mm, 
5.0  mm,  9.0  mm,  12.0  mm  and  15.0  mm.  With  the  shroud  removed,  an  infinite 
channel  spacing  was  created. 

When  measuring  the  steady  state  temperatures  of  each  heater  element,  mea¬ 
surements  were  taken  at  10  minute  intervals  and  compared  with  the  previous  data 
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run.  When  all  thermocouple  measurements  varied  less  than  0.05°C,  the  final  steady 
state  temperature  measurements  were  recorded.  The  temperature  acquisition  pro¬ 
gram  of  the  previous  study  [Appendix  A]  was  executed  to  record  and  store  data. 

Ambient  temperature  was  the  average  temperature  from  three  thermocouples 
located  in  the  fluid  near  the  top,  bottom  and  center  of  the  tank.  All  experiments 
were  stopped  if  these  three  temperatures  varied  by  more  than  0.1°C. 

When  measuring  the  fluid  temperature  with  the  movable  probe,  the  data  ac¬ 
quisition  computer  was  programmed  to  continuously  monitor  the  voltage  output 
from  the  traversing  thermocouple  probe.  This  allowed  for  point  to  point  tempera¬ 
ture  measurements  within  the  flow. 

After  all  data  was  recorded,  the  substrate  was  painted  black  to  minimize  light 
reflection.  Particles  of  pliolite  were  mixed  with  the  distilled  water  to  visualize  fluid 
flow.  The  /  stop  of  the  camera  was  set  to  2.8  and  exposure  time  for  all  photographs 
was  20  seconds.  The  camera  was  mounted  on  a  tripod,  approximately  three  inches 
from  the  glassplate  wall  of  the  tank.  Camera  height  was  adjusted  and  the  laser 
was  aligned  so  that  all  eight  heater  surfaces  would  be  in  the  field  of  view.  When 
the  shroud  was  in  position,  the  camera  single  was  inclined  slightly  to  avoid  the 
undesirable  light  reflection  off  the  shroud.  A  shroud  spacing  of  12  mm  from  the 
substrate  was  the  minimum  channel  width  at  which  useful  pictures  could  be  taken. 
Figure  2.4  shows  the  arrangement  of  the  camera,  laser  and  test  surface. 

C.  DATA  ANALYSIS 

The  convected  rate  of  energy  from  each  heater  was  determined  using  an  energy 
balance  where: 


Qconv  =  Qin  —  Qcond 


(2.1) 
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CAMERA 


Figure  2.4:  Laser  and  Camera  Arrangement  for  Photography  in  the  x-y 
Plane. 
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The  energy  in  Q,n  was  regulated  by  the  power  supply  and  determined  as 
described  in  the  initial  study.  The  energy  loss  by  conduction  was  estimated  by 
using  Fourier’s  law  of  conduction  for  one  dimension: 

AT 

Qcond  =  A,kM  (2.2) 

A ,  is  the  surface  area  of  each  heater  element,  kpg  is  the  thermal  conductivity  for 
plexiglass.  A Tavg  is  the  temperature  difference  between  the  temperature  average  of 
the  five  thermocouple  locations  on  the  heater  and  the  temperature  at  the  back  side 
of  the  substrate.  Ax  represents  the  minimum  thickness  of  the  substrate  accounting 
for  the  grooves  cut  for  wire  routing. 

A  characteristic  length,  Z,  based  on  the  geometric  dimensions  of  the  heaters 
.  was  chosen  for  determining  non-dimensional  temperature  excess,  4~,  a  flux  based 
Grashof  number,  Gr*,  and  Nusselt  number,  Nu. 

The  characteristic  length  is  the  heater  surface  area  divided  by  the  heater 
perimeter. 


I 


(2.3) 


Local  heat  transfer  relationships  were  determined  by  the  following  equations: 


"  Qconv 
q  ~  A. 


(2.4) 


(2.5) 


(2.6) 
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r.  90qL< 

Gr  =  k,*> 

(2.7) 

kj 

(2.8) 

0  =  AT 

(2.9) 

0 

NDT  =  - 
0O 

(2.10) 

The  data  reduction  program  [Appendix  A]  was  used  for  all  calculations.  Sam¬ 
ple  calculations  are  contained  in  Appendix  B.  The  various  symbols  used  above  are 
defined  in  the  List  of  Symbols. 

Fluid  properties  were  determined  from  curve  fits  listed  by  Willson  [Reference 
7].  These  are  based  on  the  data  provided  in  Incropera  and  Dewitt  [Reference  8]. 
All  properties  were  evaluated  at  a  film  temperature.  The  film  temperature  was 
determined  for  each  case  from  the  local  and  ambient  temperatures. 

A  non-dimensional  length,  was  used  in  the  presentation  of  the  data.  The 
vertical  length  £  is  the  distance  from  the  bottom  edge  of  the  substrate  to  the  local 
point  of  interest  divided  by  the  overall  length  of  the  substrate. 

Quantitive  developments  were  expressed  as  plots  of: 

•  Nu  vs.  Gr* 

•  Non-dimensional  local  temperature  excess  vs.  local  thermocouple  location. 

•  Temperature  distribution  in  the  fluid  adjacent  to  the  test  surface. 
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III.  TRANSPORT  ADJACENT  TO  AN 
ISOLATED  VERTICAL  SURFACE 


A.  FLOW  VISUALIZATION 

Photographs  of  the  steady  state  natural  convection  flow  in  the  x-y  plane  are 
shown  in  Figure  3.1.  The  laser  was  aligned  so  that  the  plane  of  light  would  pass 
through  the  geometric  center  of  the  heated  region  of  each  component.  This  figure 
presents  all  elements  at  input  power  levels  of  0.2,  1.0  and  2.0  watts  respectively  from 
left  to  right. 

Buoyancy  forces  near  the  heaters  create  a  strong  upward  flow.  The  flow  ap¬ 
pears  to  be  laminar  and  is  parallel  to  the  substrate  after  the  first  component.  A 
boundary  layer  begins  to  develop  out  from  the  leading  heater.  The  boundary  layer 
thickness  increases  to  about  2.5  component  heights  at  the  trailing  element.  At  1.0 
watt,  the  boundary  layer  is  about  2.0  component  heights  thick  and  decreases  to 
about  1.0  component  height  in  thickness  at  the  trailing  element  for  2.0  watts.  The 
fluid  velocity  increases  and  the  thickness  of  the  boundary  layer  decreases  with  in¬ 
creasing  component  power  input.  Ambient  fluid  is  entrained  below  and  along  the 
length  of  the  boundary  layer.  Figure  3.2  visualizes  the  flow  in  the  x-z  plane  at  the 
respective  power  levels  of  the  previous  figure.  The  flow  is  vertical  along  the  array 
for  about  the  width  of  the  heated  region  of  the  strips.  Weak  entrainment  from  the 
sides  is  evident.  Side  entrainment  increases  with  increasing  component  heat  input. 
Figure  3.2  clearly  shows  that  for  the  region  dong  the  array,  the  flow  can  be  treated 
to  an  excellent  approximation  as  two  dimensional. 


i 
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Figure  3.1:  Steady  Flow  in  the  x-y  Plane  for  Power  Levels  of  0.2,  1.0  and 
2.0  Watts 


B.  QUANTITATIVE 

Figure  3.3  is  a  non-dimensional  representation  of  local  temperature  excess 
levels  measured  at  all  thermocouple  locations  for  an  input  power  level  of  2.0  watts. 
As  expected,  the  measured  component  center  temperatures  were  the  highest.  The 
top  edge  temperatures  always  exceeded  the  temperatures  at  the  other  edges.  The 
left  and  right  edge  temperatures  differed  by  a  maximum  of  22%  at  2.0  watts.  The 
heater  strip  had  an  11  mm  tab  on  the  left  side  for  power  lead  connections.  Left  edge 
conduction  losses  were  therefore  expected  to  be  greater  than  the  right  side,  resulting 
in  slightly  lower  temperatures  on  the  left.  After  the  third  component,  the  bottom 
edge  temperatures  exceeded  the  excess  temperature  of  the  sides.  The  side  edges 
benefited  more  from  fluid  entrainment  than  the  bottoms  which  are  in  the  wake  of 
the  heated  upstream  components.  In  all  experiments,  the  right  edge  temperature 
for  component  seven  was  lower  than  the  established  trend.  It  is  believed  that  the 
thermocouple  placement  at  this  location  may  have  slipped  when  the  heater  foil  was 
pressed  into  place. 

Heat  transfer  coefficients  were  computed  from  the  data  as  outlined  in  equation 
2.6.  They  indicated  a  downstream  decreasing  trend  at  a  given  power  level.  This 
is  expected  since  the  non-dimensional  temperature  excess  is  inversely  proportional 
to  h.  As  seen  in  Figure  3.3,  the  general  trend  tends  to  be  increasing  downstream 
excess  temperatures. 

As  the  component  heat  flux  is  increased,  the  buoyant  flow  becomes  more  vig- 

\ 

orous.  The  component  temperature  levels,  as  a  result,  decrease.  The  calculated  heat 

\ 

transfer  coefficients  are  greater  than  those  for  the  same  component  at  lower  power 
levels.  With  the  larger  values  of  h  at  higher  power  levels,  there  is  a  corresponding 
decrease  in  the  non-dimensional  temperature  excess. 
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Figure  3.3:  Non-dimensional  Temperature  Excess  vs.  Position:  2.0  Watts 
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Figure  3.4  is  a  logarithmic  scatter  plot  of  Nu  vs.  Gr*.  The  cluster  of  data 
points  at  approximately  a  single  Gr *  represents  the  Nusselt  numbers  for  compo¬ 
nent  centers  at  each  power  level.  The  Nusselt  numbers  decrease  for  downstream 
components,  due  to  the  increasing  component  temperatures. 

As  in  the  previous  study,  a  modified  Nusselt  number,  N u*,  was  calculated  by 
using  a  scaling  factor  to  account  for  the  increase  in  convected  energy  of  the  flow  at 
downstream  locations. 


Nu*  =  Nun  x  S,!4  (3.1) 

Nun  and  Sn  are  the  local  Nusselt  numbers  with  the  respective  heater  element  num¬ 
bers. 

This  modified  data  is  plotted  in  Figure  3.5.  The  linear  least  squares  fit  through 
the  data  with  a  ±5%  band  is  also  presented.  Using  the  scaling  factor  above,  nearly 
all  of  the  modified  data  falls  within  this  band.  Based  on  this  curve  fit,  the  following 
heat  transfer  relationship  for  the  given  array  geometry  was  obtained: 


Nu *  =  1.115(Gr*)' 


for  250  <  Gr*  <  7000  and  Pr  ~  7. 


C.  EFFECT  OF  COMPONENT  SPACING  ON  HEAT  TRANSFER 

Now  that  a  heat  transfer  correlation  has  been  developed  for  the  given  geometry 
of  the  array,  it  is  desired  to  see  how  a  change  in  component  spacing  would  affect 
these  relationships.  This  was  achieved  by  powering  every  other  and  every  third 
heater.  This  had  the  same  effect  as  heater  elements  spaced  50.8  mm  and  76.2  mm 
center  to  center.  The  resulting  variations  of  non-dimensional  temperature  excess  at 
the  heater  centers  are  shown  in  Figures  3.6  and  3.7. 
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Figure  3.4:  Data  for  Nusselt  Number  V3.  Flux-Based  Grashof  Number 
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Figure  3.5:  Modified  Nusselt  Number  vs.  Flux-Based  Grashof  Number 
with  Least  Squares  Fit 
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Figure  3.6:  Non-dimensional  Temperature  Excess  vs.  Position  with  Four 
Elements  Powered 


Figure  3.7:  Non-dimensional  Temperature  Excess  vs.  Position  with 
Three  Elements  Powered 
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When  heaters  2,  4,  6  and  8  are  powered,  there  is  a  slight  decrease  in  the  non- 
dimensional  temperature  at  the  first  downstream  component  beyond  the  leading 
heater.  Further  up,  component  temperatures  vary  little  from  that  of  the  second 
heater.  When  heaters  2,  5,  and  8  are  powered,  the  first  two  component  temperatures 
remain  unchanged.  The  trailing  heater  has  a  small  drop  in  its  non-dimensional 
temperature. 

For  the  increased  component  spacings  in  Figures  3.6  and  3.7,  the  convected 
energy  due  to  upstream  heaters  has  only  a  weak  effect  on  downstream  components. 
The  fluid  momentum  resulting  from  the  buoyancy  forces  actually  provides  moderate 
improvement  in  heat  transfer  for  the  downstream  heaters.  Heat  transfer  relation¬ 
ships  were  obtained  for  all  components  by  treating  the  elements  in  the  array  as 
single  heaters.  The  resulting  component  based  Nusselt  numbers  became  almost  in¬ 
dependent  of  the  downstream  component  position  in  the  expanded  array.  These  are 
plotted  against  the  flux-based  Grashof  number  in  Figures  3.8  and  3.9.  The  least 
square  fit  to  the  correlated  data  for  all  unshrouded  heaters  is  shown  for  comparison 
in  each  figure. 

D.  TEMPERATURE  MEASUREMENTS  WITHIN  THE  FLUID 

Measurements  of  the  thermal  transport  in  the  fluid  adjacent  to  the  test  surface 
were  made  using  the  movable  thermocouple  probe.  The  thermocouple  junction 
was  initially  centered  on  each  component  and  then  moved  outward  in  the  direction 
normal  to  the  substrate.  Temperatures  axe  seen  in  Figure  3.10  for  input  power  levels 
of  0.2  and  1.0  watts.  The  thermal  boundary  layer  in  both  cases  is  well  within  the 
momentum  exchange  region  which  is  approximately  1.6  cm  thick  for  component  8 
at  1.0  watt.  This  is  expected  due  to  the  relatively  large  Pr  of  about  7.0  for  water. 
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Figure  3.8:  Nusselt  Number  vs.  Flux-based  Grashof  Number  with  Four 
Elements  Powered 
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Figure  3.9:  Nusselt  Number  vs.  Flux-based  Grashof  Number  with  Three 
Elements  Powered 


Figure  3.10  suggests  that  it  may  be  possible  to  move  a  shroud  wall  to  within  one 
component  height  without  changing  the  rate  of  heat  transfer  significantly. 

The  traversing  probe  was  slowly  swept  spanwise  in  the  z-direction  past  the 
centers  of  the  heaters  3,  5  and  8  at  a  distance  of  1.0  mm  outward  in  the  y-direction 
from  the  heater  surface.  Resulting  temperature  variations  are  shown  in  Figure  3.11 
for  a  power  level  of  1.0  watt.  The  component  center  is  referred  to  as  0.0.  A  scale 
line  provides  a  reference  for  the  width  of  the  heated  region  of  each  component.  A 
large  drop  in  temperature  is  seen  to  occur  beyond  Z  =  ±6  mm.  For  the  center  50% 
of  each  heater,  there  is  relatively  small  change  in  temperature,  again  suggesting  a 
largely  two  dimensional  flow  there. 


o 

§ 20 
WW 
o  cj 

05 

D 

CO 


W  W 
CJO 
<5  ■< 
W  W 
05  05 
DD 
COCO 


WWW 

CJOCJ 

<<< 

WWW 

05  05  05 
DOO 
COCO  CO 


W  Wl 
CJ  CJ  i 

<  < 

W  Wl 

05  051 

dd: 

COCO  i 


□  o  <  +  x  o  o  a  x 


li  *3 


i  i  i  i  i  i  i  i  i 

i  i  i  i  i  i  i  i  i 

■  -i - 1- -i - 1- -i - - i.B-i- 

. . 

i  i  i  i  i  i  K  i  ® ' . .  -  -H 

•  - ' - 1__( - i  -  a-  -  gi, 

i  i  i  iSS  i  BT  u  -  •®  . . .  <;  ®V 

._i - -gj®jg  4>j  -.V.'Cj&jp-X. 

1  •  > — r '  '  -L- '  i  Ji . 

i  i  pi  tx©  i  i  •jB  i  m,i 

•  -  >-  -  L  -a  — 1  7-/ 

^tlsw'  ff''  1/  (t 


0  1  SB  0B  S'!  0'E  SI  01  SO  00 

(o  oaa)  aaaxvaadKax  vxaaa 


Figure  3.10:  Fluid  Temperature  Distribution  in  the  y-Direction  for  0.2 
and  1.0  Watts 
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IV.  TRANSPORT  IN  A  VERTICAL  CHANNEL 


A  plexiglass  shroud,  simulating  the  back  side  of  a  neighboring  circuit  board, 
was  positioned  to  form  a  vertical  channel  of  various  widths. 

A.  FLOW  VISUALIZATION 

Useful  photographs  could  only  be  taken  with  the  shroud  positioned  no  closer 
than  12  mm  from  the  substrate.  All  visualizations  presented  are  at  power  levels  of 
0.2,  1.0  and  2.0  watts  respectively.  Figure  4.1  shows  the  natural  convection  flow  in 
the  x-y  plane  with  the  shroud  positioned  12  mm  from  the  substrate.  The  shroud 
only  allows  entrainment  from  the  sides  and  bottom.  Initially  the  flow  is  not  parallel 
to  the  channel  walls.  The  entrained  fluid  entering  from  the  bottom  has  velocity 
components  in  the  x  and  y  directions.  As  the  boundary  layer  develops  to  the  width 
of  the  channel,  the  flow  is  straightened  and  conforms  to  the  channel.  This  same  flow 
behavior  occurs  when  the  channel  width  is  increased  to  15  mm  as  shown  in  Figure 
4.2  The  effects  of  entrainment  are  more  noticeable  downstream  since  the  boundary 
layer  has  more  time  to  develop  with  the  wider  channel. 

Figure  4.3  presents  flow  visualization  in  the  x-y  plane  with  the  shroud  posi¬ 
tioned  6  mm  from  the  substrate.  These  photographs  were  taken  with  the  camera 
aligned  in  the  direction  normal  to  the  surface.  With  the  shroud  in  position,  side 
entrainment  is  much  more  pronounced  than  the  case  without  a  shroud  [Figure  3.2). 

B.  QUANTITATIVE 

Figures  4.4  and  4.5  are  non-dimensional  representations  of  temperature  excess 
levels  for  all  center  located  thermocouples  for  each  shroud  condition.  These  two 
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Figure  4.1:  Steady  Flow  in  the  x-y  Plane  for  Power  Levels  of  0.2,  1.0  and 
2.0  Watts  with  12  mm  Shroud 


Figure  4.2:  Steady  Flow  in  the  x-y  Plane  for  Power  Levels  of  0.2,  1.0,  2.0 
Watts  with  15  mm  Shroud 
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figures  depict  the  spread  in  data  on  the  same  plotting  scale  for  the  maximum  and 
minimum  power  levels. 

Again  it  is  noted  that  the  general  trend,  regardless  of  input  power  level,  is 
increasing  downstream  temperature  excess.  As  power  is  increased  from  0.2  watts 
to  2.0  watts,  the  values  for  the  non-dimensional  temperature  excess  decrease.  As 
expected,  the  heat  transfer  coefficient  increases  for  corresponding  components  as 
the  power  level  increases. 

There  is  little  difference  in  the  non-dimensional  temperature  excess  as  the 
shroud  is  positioned  6  mm  or  more  from  the  substrate.  When  the  shroud  is  closer 
to  the  substrate,  it  reduces  the  ambient  fluid  entrainment,  thereby  increasing  com¬ 
ponent  temperatures.  This  is  less  pronounced  at  the  higher  power  levels  as  seen  in 
Figure  4.5. 

Figure  4.6  shows  the  actual  temperature  difference  for  the  component  centers 
versus  the  shroud  position  at  2.0  watts.  For  this  case  and  ail  other  power  levels, 
the  actual  temperature  difference  was  greater  when  the  shroud  was  closer  than  6 
mm  from  the  substrate.  When  positioned  6  mm  or  more  from  the  test  surface,  the 
actual  temperature  difference  was  within  0.5°C  of  the  unshrouded  levels. 

The  modified  Nusselt  numbers  for  component  centers  were  calculated  using 
the  same  correlation  as  in  the  previous  chapter.  These  values  were  compared  on  a 
logarithmic  plot  with  the  flux-based  Grashof  numbers  at  each  shroud  position.  The 
base  line  curve  for  the  linear  least  squares  fit  of  the  unshrouded  data  with  ±5% 
scatterband  was  also  plotted  for  comparison.  The  scatter  plot  for  data  when  the 
shroud  was  positioned  1.5  mm  from  the  substrate  is  shown  in  Figure  4.7. 

The  only  points  in  close  agreement  with  base  line  data  are  those  points  corre¬ 
sponding  to  the  leading  heater  element.  Regardless  of  power  level  or  shroud  posi¬ 
tion,  the  leading  heater  can  be  treated  as  a  single  discrete  heat  source  independent 
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Figure  4.4:  Non-dimensional  Temperature  Excess  vs.  Position  with 
Shroud  at  0.2  Watts 
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Figure  4.5:  Non-dimensional  Temperature  Excess  vs.  Position  with 
Shroud  at  2.0  Watts 
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Figure  4.6:  Temperature  Difference  vs.  Shroud  Spacing:  2.0  Watts 
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of  downstream  array  geometry.  A  scatter  plot  of  data  for  the  shroud  position  of  3.0 
mm  is  not  much  better  than  the  previous  case. 

Figure  4.8  shows  a  plot  of  modified  data  when  the  shroud  spacing  is  increased 
to  6.0  mm.  The  same  trends  seen  here  were  found  at  all  larger  spacings.  All  data 
are  seen  to  be  in  close  agreement  with  the  correlating  equation  thus,  with  a  channel 
width  of  6.0  mm  or  more  in  width  and  for  the  given  geometric  array,  the  modified 
Nusselt  number  can  be  calculated  to  within  5%  accuracy  using  Equation  3.2  for 
250  <  Grl  <  7000  and  Pr  ~  7. 
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Figure  4.8:  Modified  Nusselt  Number  vs.  Flux-based  Grashof  Number 
with  Shroud  at  6.0  mm 
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V.  CONCLUSIONS 


Flow  visualization  and  quantitative  analysis  provided  several  significant  con¬ 
clusions: 

•  The  flow  was  laminar  over  the  given  range  of  power  levels. 

•  The  convective  flow  is  nearly  two-dimensional  when  the  shroud  is  not  present. 
This  condition  greatly  simplifies  computational  modelling  of  such  flows. 

•  A  thermal  layer  was  embedded  within  the  momentum  region.  This  la>er  be¬ 
came  thinner  at  increasing  levels  of  heat  flux. 

•  A  general  trend  of  increasing  component  temperatures  downstream  at  each 
power  level  was  found. 

•  No  significant  gain  in  the  rate  of  heat  transfer  for  channel  widths  of  6.0  mm 
or  more  was  found  over  the  unshrouded  configuration. 

•  For  geometrically  similar  arrays,  heat  transfer  relationships  for  downstream 
components  can  be  determined  by  using  the  correlating  equation  for  this  ge¬ 
ometry. 

•  As  component  spacing  is  increased  to  two  times  or  more  of  that  for  the 
present  configuration,  upstream  effects  become  negligible  and  the  heaters  can 
be  treated  as  an  individual  discrete  heat  source. 
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VI.  RECOMMENDATIONS 

In  continuation  of  this  study,  it  is  suggested  to  use  a  similar  array,  however: 

•  Study  the  effects  of  even  smaller  component  spacings  on  heat  transfer. 

•  Flank  each  side  of  the  single  array  with  a  column  of  flush  mounted  heat  sources. 

•  Spray  the  heater  elements  with  temperature  sensitive  liquid  crystals. 

•  Measure  the  variation  of  temperature  within  the  thermal  layer  at  higher  power 
levels  for  all  components  and  spaces  between  components. 

•  Explore  the  transient  evolution  of  transport. 


APPENDIX  A 
SOFTWARE 

A.  TEMPERATURE  ACQUISITION  PROGRAM 
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5u 

REAL  Voltsi S3  i 

60 

REAL  Tempi  59/ 

70 

CREATE  8DAT  “ 

•  ,700  ,0  ,  ‘  ,48 

60 

ASSIGN  9Patm3  TO  “ 

: .700 .0.  ' 

50 

PRINT  ” 

SURFACE  #6" 

100 

PRINT 

1  10 

OUTPUT  70S; 'CGNrMEAS 

DCU.100.USE  0" 

120 

ENTER  709;Voltsio0 ) 

1  33 

OUTPUT  70S ; " CONFMEAS 

OCU, 100- 105, USE  0" 

140 

FOR  1-0  TO  5 

150 

ENTER  709 (Voltsi  I  > 

l  S3 

T em pi  I  /-. 0006797+125325. 1 328*Vci t  s 1  I  1  >-(607739. 2467*1  Volts;  I  >  *Vo 1 1  s (  I  >  >  )- 

21952 

034. 3364»i  Uolt  si  I  >  *‘  3  )  >+<  8370810996. 1  974*1  Voltsi  I>''4  )  >+.075 

170 

180 

PRINT  “T.  C.  #“;I+t 
NEXT  I 

Volts  O.C.  ”  ; Uolt si  I > ."Temp.  OEo.  C  "i Tempi  I) 

193 

PRINT 

230 

PRINT  " 

SURFACE  47" 

210 

PRINT 

220  OUTPUT  709;  "CONFNE.-r  DCV  ,  1 06-1 1  I  .USE  3  ’ 

233  Fur.  1-5  TO  I  ! 

0*13  ENTER  739:VoI;s(I/ 

253  Tempi I  )• . 0006797+; 25325. 1 328* Vo I t s ( I > )-( 537789.2457 ♦iUoits(  I  )*Vol t  s( I  ) )  )-< 
21 952 034. 3364 •( Volt  si I  )"3 i )+< 33708 10996 .1874*1  Vo itsil>'4>  >+.075 


260 

PRINT  “T.  C.  *";I+1 ,* 

Uolts  O.C.  " ; Vo 1 1 s i I ) , "Temp .  DE6.  C  ";Temp(I> 

270 

NEXT  I 

2S3 

PRINT 

290 

PRINT  - 

SURFACE  #6" 

300 

PRINT 

310 

OUTPUT  709;  'CONFllEAS  OCU 

,l  12-1 17, USE  0" 

320 

FOR  I • 1 2  TO  17 

333 

ENTER  709  ;  Uo  J.  t  s  i  I  > 

340 

Tempi  I  )-.0006797+i  25825. 

1 328*Voi t si  I >  >-( 807799 .2467*<Volts< I )*Vol t s(  I > / > 

2 1  552034 . 3364*1 Vol t si  I >"3 > >+<  83708!  0S9S  .  1874*t Voltsi I )*4>  )+.07S 

353  PRINT  "T.  C.  t" ;  1  +  1  ,  *  Volts  O.C.  “  i  Vol  t  s<  I  > ,  “Temp  .  OES.  C  "sTempil) 

35C  NEXT  I 
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I 

! 


CF.INT 

35C 

rr.  Ii’vT 

35C 

FF.KJT 

*GC 

jL1  :  r  J  :  .  C  r  •  ‘ir 

*♦  •  G 

:  •%  t  *  :  5  ■  9 

*10 

ENTER  70=  :  ,'Cit  s  I  . 

430 

Temp; 1  )• .  0005757+' 

2 1  952034 .3564* (  Vo  I  t s(  I  ) " 

440 

PRINT  “T .  C.  #";I 

450 

NEXT  I 

460 

OUTPUT  709; -CONFME 

470 

FOR  1-20  TO  23 

4S0 

ENTER  703 t  Vc 1 1  s (  I  ) 

4  90 

Temp* I  /-.G00E737-' 

Z :  9  91 

C“4 . ZZE 4 * • , zi ’ ;  I 

500 

PRINT  “ T  .  C  .  s ;  I 

5 1 0 

NEXT  I 

520 

PRINT 

530 

PRINT  " 

540 

Pn  I  NT 

E9C 

OyTF'j?  “Cs ; ''  IDNF'-'E 

550 

:nc  ; *24  T0  ZS 

ENTER  705;  volt  S'.  I  . 

5S0 

Temp  •  I . 0005 “57*i 

21953 

034 . 3354*  t  Uol  is*.  I  > 

550 

PRINT  "7.  C.  S  ;  I 

500 

NEXT  I 

9  i  0 

PRINT 

520 

PRINT  " 

£30 

=R:I  NT 

540 

OUTPUT  70S : " CONFNE 

B5G 

FOR  1-30  TC  35 

55C 

ENTER  “03; vc it  s'  I  • 

570 

7  *  rr  ~  I  *  .  O00E7B71*" 

Z  ;  39Z 

£“>4.3354*  Uoits-I1 

□  50 

PRINT  'T.  C. 

690 

NEXT  I 

700 

FOR  J-l  TO  14 

7i0 

PRINT 

Uolts  D.C.  ‘ ; Uol t s ( I  )  ,  “Temp .  DEG.  C  “;Temp(I) 


;  :  emp 


SUFFICE  S4" 


3  .  JZZ.  i> 


537061  055E .  i  £74 *  i  Vo  1 1  s  i  I  / '  4  •+  .  C“5 
, "  U  c  i  t  s  D.  3 .  "  5 volts'  I  >  ,”Temp  .  DEG.  C  "sTemo1.  I; 


SUFFICE  S3" 


£3 TOE i SS35 . 1  £74* . vo i t  s  . 


IE~*  '-’c Its'  I  <«Uslt8 • 
4 • it . GTE 
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720  NEXT  J 

730  PRINT  "  SURFACE  42" 

740  PRINT 

7S0  GUTPUT  709;  “CONFMEAS  DCV  ,2 1 S-Z 1 9  .USE  0" 

750  FOR  1-35  TO  39 
770  ENTER  70S i Volt I) 

7  50  Tefxp'  I  .•».  00C'5’'S7 *• 25325 . 1 3ZE -Vc i t  s . I  1  •  -  .  5G775  9 .  2457*  <  Vc  i  t  s  •  i  ■  ♦  vc  1 1  s  I  '  1  •- 
3 ‘ 952034 , 3354  *  <  Volts  I1  3  ■  •' 32705 ' 293c  .  1  5  ~4  » •.  vo  i  1 1  ■  I  ?  4  '  +  .  C” 

750  PRINT  T.  C.  *“  s  1*1  . "  Volts  D.C.  "  iVrlts  I  ' ,  '7e"ip.  DEE.  2  '  ;7ewp-.  I  .■ 

SCO  NEXT  I 

610  OUTPUT  70S; "CONFMEAS  DCV  ,300-30 ’  .USE  0" 

S20  FOR  1*40  TO  41 
S30  ENTER  709; Volts; I  ) 

840  Temp(  I >-.00067974(25825. 1328*Volts( I >  >-< 507789 . 24E7*( Vol t s< I >*Volts; !)>>-< 
21952034.3364*' Volts! I >*3 > >+( 93708 1 0995 . 1 874* i Volts; I >"4  >  >-.075 


350 

PRINT  “T.  C.  #" ;I+1  , 

Volts  D.C.  “ ; Vol t s( I  /  .  “Temp ,  DE5.  C  " ; T  emp ( I 

550 

NEXT  I 

370 

PRINT 

550 

PRINT  ” 

SURFACE  si" 

550 

PRINT 

900 

Output  70S; "COnfmeas 

DC V, 302 -307  . USE  0" 

<m  .  <s 

3  I  <C 

FOR  1-42  TO  47 

920 

ENTER  709 ; Vol t  a ; I > 

930 

Temp' I )- , 0006 7 97+ ( 255 

25. 1 328* Vol t s  I >  507755 . 2457* (Volts; I >»Volts; I 

2!  952 

034.3364*;  Volts;  I  >'  3  > )+'.  8370510996. 1374* 'Volt  S'.  I  >*4  >  >+.075 

540 

FFINT  "T,  C.  t" ; 1+1 , 

Volts  D.C.  "; Veits! I >, "Temp.  DEE .  C  "(Temp'I 

350 

next  : 

550 

'OUTPUT  @RatnZ;Temp; * 

570 

p=:nt 

9=0 

PRINT 

930 

RRINT  " 

BATH  TEMPERATURES  \ TOP  TO  BOTTOM  1" 

1  00G 

FPJNT 

1010 

OUTFUT  709; "CONFMEAS 

DCV. 3)7-31 9, USE  0" 

1020 

FOR  1-57  TO  59 

1  030 

ENTER  709 ; Vc its;  I  > 

1  £40 

Temp (  I }— .00G5797+ (255 

25 .  1 32 a* Vo  1 1  s '  I  > 50-7SS .  2467* !  Vc  1 1  s '.  I  )* Vo i  t  s  '  I 

2 1 952034 . 3354*<  Vci t  s ; 1 >  '3  > >+( B3703 1 0995 . l874*(Volts; I  )  4  >  »+ . 075 
1050  PRINT  “VOLTS  D.C.  "; Volts; I >  ."TEMP.  DEG.  C  ";Temp(I> 

1050  NEXT  I 
1 070  ENO 
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B.  DATA  REDUCTION  PROGRAM 


i  2  I  i  I  •  I  ;  .  •  '  i  :  :  !  i  i  I  i  i  .  :  .  i  ■  • . :  \  ‘  i 

20  1  ! 

32  i*  DATA  REDUCING  PROGRAM 

43  t  i 

S3  II  |  I  M  ;  II  i  I  I  I  i  M  I  I  I  i  !  I  I  I  ;  I  !  M  I  I  I  i  i  l|  I  I  l  M  I  I  M  I  I  I  I  I  :  II  I  !  U  M  I  i  I  i  I  I 

52  REAL  T  {  45  ■  ,  Ocondt  S  ' 

70  PERL  Power  .Tamo  ,S  ,Kpa  ,0ei.<  , As  .Per  ,Ta  .Qccnv  ,Lbar 

30  REAL  Qf lux  ,Deit  .Tfilm .Beta ,5pvol .Dynvisc 

90  RERL  Kinvisc  ,Kf  ,Pr  .Theta  ,Ndt  ,H>.  .Grmcd  ,Lnu  ,Lgr 
100  PRINTER  IS  1 

1  10  PRINT  "INPUT  DRTR  FILE  NAME" 

122  INPUT  Filenames 

130  PRINT 

140  PRINT  "INPUT  POWER  VALUE  IN  WATTS” 

150  INPUT  Power 

150  FRINT 

170  FRINT  "INPUT  SHROUD  SPACING  IN  mm" 

180  INPUT  Spaces 

190  FRINT 

200  PRINT  "INPUT  AMBIENT  TEMPERATURE" 

210  iNFUT  TamQ 

220  FFInTER  IS  ~3  i 

230  PRINT  "POWER  LEUE-  \ WATTS)  •" ,Pcwer 

252  PRINT  "SHROUD  POSITION  ■  nn ;  Spaces 

260  PRINT 

270  PRINT  "AMBIENT  TEMPERATURE  >' C  >  :”,Tamc 
280  PRINT 

290  PRINT 

300  ASSIGN  ®Path  TO  Filenames 

310  ENTER  9Pat hi T( *  ) 

323  G-9.81 

330  Kpg* .142! 

340  Del.*.* .  00673 1 

350  As= . 0335* . 3C75 

363  rer»2* .3239+2 • .0075 
370  Lbar»As/ Per 

380  A$*"Qcond" 


;  | 

i  i 
i  i 
i  i 


390  B$» "  >.  'l  " 

4-30  Cs-'NO"' 

410  D$» ''LOS'  Nu •.  ■  ” 

420  ES»"L0G\Grv*  >" 

430  ri«"DELT“ 

440  GS-'H,.  “ 

450  HS- " Nux " 

460  I$-"Grx«“ 

470  | j {{!{!!!!{!! I !!!<!{!!!!!! I  i  I  M  I  I  {  f  |  i  II  |M  m  m  i  I 

480  !  COMPUTE  CONDUCTION  LOSSES  ! 

4  90  M  I  |  M  |  m  )  )  |  i  I  I  |  !  i  !  n  !  M  i  !  1  1  1  >  i  i  !  !  !  1  1  !  !  M  !  !  I  i  M  i 

500  J-8 

510  FOR  1-0  TO  42  STEP  6 

520  T  a  -  (  T  <  I  >  +  T  :  I  + 1  >+T\I+2>+TfI+3>+T\I+4,)/5 
530  uccnoi  j  >-As«Kpg*i  Ta-T;  1*5  •  ■,  D»i> 

540  J-J-l 
550  NEXT  I 

550  i  m  I  M  !  !  I  !  M  1  1  M  1  !  I  i  M  !  M  *  i  i  M  1  !  *  !  !  !  M  !  M  !  !  i  M  !  I 

570  N-0 

530  FOR  J-8  TO  I  STEP  -1 
550  Cm  =0 

60S  J-!  '*25 . 4t?4  •  3 '  T 

510  X  o  - 1  J  - 1  )  *  2  5 . 4  *  6  6  >  •  3  ’  7 

620  Xc-C C  J- 1  >*25.4*70 »/3i 7 
630  PRINT 
640  PRINT 

650  PRINT  ''  SURFACE  #"  ,J 

660  PRINT 

670  PRINT  USING  "2X  .AAAAA ,4X ,AAA ,4X ,AAA ,4X .AAAAAAAA  ,2X .AAAAAAAAA ,3X , AAAA  ,5X  ,AA 
,6X  , AAA  ,4X  , AAAA “ ;A$  ,BS  ,CS  ,DS ,ES ,FS ,GS  ,H$ , I* 

630  FOR  I-N  TO  N+4 
630  Qconv-Pouei — Gcond( J  ) 

700  Of  lu.- -Oconv/As 
7 1 0  Del t -T( I  i-Tamb 
720  Tfiim-273. 1 5+i T. I  >+Tanb  >/2 

730  6eta«\ -7 . 944856E-3*<  Tf i in" 2  >  >+< 5 .735S479E-5*Tf i In  >-.0097310563 

740  Spvol-<  4.69699E-9*< Tfiln'2 ) >+( -Z .53745E-6+Tf Un >+.001341903 
750  Dynvi  sc-(  3.25485 1  1  E-7*<  Tf  i  ln"2  >  >+<-. 0002 1 4744S7«Tfi In >+. 036 1 66792 

760  Kinvisc-Spvol*Oynvisc 
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770  Kf-< I .41318181 S£-3»Tfilm >* . !3oS 

780  Pr-( 4 . B5706E-3* < Tf i lm  2  -  >  +  i -2 . 322084*1 f ilm >+4E3 .33 1 9 
790  Theta«Qf lux*Lbar/Kf 
800  Ndt-Dslt/Theta 
010  Hx-Of lux/Dei t 

620  Srmod-( G*Beta*Qflux*(  Lbar  )'4  )/( Kf  *Kmvisc''2  ) 

830  Nux“H>  ♦Lbar./K  f 
840  Lflr«LoTt 6rnod ) 

850  Lnu”L3T< Nux ) 

8B0  IF  Cnt“0  THEN  Xtol“Xt 

870  IF  Cnt-1  THEN  Xtol-Xb 

880  IF  Cnt > 1  THEN  Xtol-Xc 

390  Cnt“Cnt+l 

900  PRINT  USING  " 1 X  ,D.DDDD  ,2V  ,D . 000 ,2X  .D.DDD  ,4X  ,0 . ODD  ,7X  ,D .COD ,4> ,00 . DD0 
DO  .0  ,2X  .OD.DDD  ,2X  ,DDDC.D“  sOccnd-' J  :  .Xtci  ,Ndt  ,Lnu  .Lgr  Delt  H>  Nu*  5-mod 

SI0  NEXT  I 


=20 

N»w+ 

530 

NEXT 

540 

END 

APPENDIX  B 
EQUATIONS 

A.  SAMPLE  CALCULATIONS 

These  calculations  axe  baaed  on  temperatures  from  data  corresponding  to  sur¬ 
face  eight  with  adl  elements  powered  to  1.0  Watt  and  no  shroud. 

1.  CHARACTERISTIC  DIMENSIONS 

P  =  2(0.0078  +  0.0239) 


0.063m 

(B.l) 

(0.0078)(0.0239) 

1.864  x  10“4m2 

(B.2) 

1.864  x  10-4 

0.0634 

2.940  x  10*3m 

(B.3) 

2.  CONDUCTION  LOSS 


T 

A  aVg 

30.40  +  27.98  +  27.09  +  27.39  +  31.11 

5 

=  28.79°C 

(B.4) 

Qcond 

(1.864  x  10“4)(0.1421)(28.79  -  26.09) 
.006731 

=  0.0106W 

(B.5) 

3.  CONVECTED  HEAT  FLUX 


Qconv  =  1.0  —  0.0106 
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=  0.9894 VT 
0.9894 
1.864  x  10-“ 


4.  FLUID  PROPERTIES 


m  31.11  +20.55 

Tfum  =  - - - h  273.15 

=  298.98# 


0  =  (-7.945  x  10~®)(298.98)  +  (5.736  x  10“5) 
(298.98)  -  9.781  x  10“3 


=  2.665  x  10-“/r 


t;  =  (4.697  x  10"9)(298.98)2  -  (2.537  x  10~6)(298.98) 
+1.342  x  10'3 

772^ 

=  1.003  x  10‘3— 
kg 

H  =  (3.235  x  10_7)(298.98)2  -  (2.147  x  10~4)(298.98) 


+0.0362 


,4Ns 


=  9.264  x  10"4  — 
rrr 

v  =  (1.003  x  10"3)(9.264  x  10"4) 

2 

=  9.262  x  10-7 — 
s 

kf  =  (1.418  x  10-3)(298.98)  +  0.1866 
W 

=  °-611^ 

5.  NON-DIMENSIONAL  TEMPERATURE  EXCESS 


9  =  31.11  -20.55  =  10.566C 


(B.6) 


(B.7) 


(B.8) 


(B.9) 


(B-10) 


(B.ll) 


(B-12) 


(B.13) 


(B.14) 
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00 

NDT 


(5308)(2.940  x  10“3) 
0.611 


10.56 

25.54 


=  0.413 


25.54°C 


6.  HEAT  TRANSFER  COEFFICIENT 


(B.15) 

(B.16) 


hx 


5308  W 

31.11  -  20.55  ~  502‘  m2K 


7.  FLUX-BASED  GRASHOF  NUMBER 


(B.17) 


,  (9.81)(2.665  x  10-4)(5308)(2.940  x  10~3)4 

Grx  ~  (0.611)(9.292  x  10-7)2  _  1965 

8.  NUSSELT  NUMBER 


(B-18) 


(502.7)(2.940  x  10‘3) 
U~  0.611 


2.419 


(B.19) 
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